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Q . Abstract 

f^ ' The recent detection of 7-rays with energy up to 10 TeV from 

^O . dense regions surrounding some Supernova Remnants (SNR) provides 

^^ I strong, though still not conclusive, evidence that the nucleonic compo- 

i-^ ' nent of galactic Cosmic Rays is accelerated in the supernova outflows. 

I ' Neutrino telescopes could further support the validity of such scenario 

2 . by detecting neutrinos coming from the same regions. We re-evaluate 

'^ I the TeV range neutrino-photon flux ratio to be expected from pion 

C^ ■ decay, finding small differences respect to previous derivations. We 

^ . apply our results and the recent HESS measurements of the very high 

K/i I energy 7-ray flux from the molecular cloud complex in the Galactic 

^ • Centre, to estimate the neutrino flux from that region discussing the 
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detectability perspectives for Mediterranean Neutrino Telescopes. We 
also discuss under which conditions neutron decay may give rise to a 
significant TeV antineutrino flux from a SNR embedded in Molecular 
Cloud complexes. 

1 Introduction 

The existence of a correlation between the position of supernova remnants 
(SNRs) and that of giant molecular clouds (GMCs) is supported by theoret- 
ical arguments and observational evidences. Infrared observations confirmed 



1 



the theoretical expectation that stars, especially massive ones, form mainly 
in the womb of GMs. Having a relatively short timelife, massive stars can 
complete their evolution and explode within the parent GMC before this is 
swept-out by stellar winds. SN explosions within, or in the nearby of, GMCs 
may trigger the formation of new stars explaining the bursting star formation 
activity observed in several dense regions of the Galaxy, e.g. in the Orion 
giant cloud, and in starburst galaxies. 

From the point of view of high energy astrophysics, SNRs interacting 
with MCs are very interesting systems. Since SN explosions are expected 
to be the main source of the bulk of cosmic rays (CRs) in the InterStellar 
Medium (ISM) jjj, a MC which is embedding, or just lay in the nearby of, 
an active SNR, should be crossed by a flux of relativistic particles which 
is signiflcantly larger than in the rest of the Galaxy. Furthermore, due to 
their relatively high gas density {riH ~ 10^ -7- 10^ cm~^), GMCs should act 
as efficient astrophysical beam-dumps converting a relevant fraction of the 
primary particle energy into secondary 7-ray and neutrino emissions (see e.g. 
|21 El)- In the following we will generally call SNRs interacting with MCs 
as embedded SNRs, also in those cases in which the SNR is not completely 
surrounded by the molecular gas. 

Evidences of a correlation between 7-ray emission and dense molecular 
gas have been found by several satellites and particularly by EGRET (see |lj) 
and have been extensively studied by many authors (see e.g. |3]). Those ob- 
servations, however, may generally be explained without invoking any local 
excess in the CR energy density respect to the value observed on the Earth. 
Only recently, a new generation of atmospheric Cherenkov telescopes, espe- 
cially HESS, were able to find the smoking-gun of the interaction of CRs 
accelerated by an active SNR with a close MC. The most interesting HESS 
observations in this respect are those of RXJ1713. 7-3946 and the very 
recent measurements of the Galactic Centre (GC) ridge emission [Zj. 

Several arguments suggest that the origin of the very high energy (HE) 
radiation coming from these sources may be hadronic. One of the most solid 
among these arguments is based on the observation of strong (up to several 
milliCauss) magnetic fields in MCs [S] (see [Hj for the GC region). Such 
strong field should give rise to dramatic electron synchrotron energy losses 
above the TeV suppressing a possible HE Inverse Compton contribution to 
the total 7-ray emission from those regions (see e.g. [Zj). Although the 
previous argument is rather convincing, a more direct evidence that HE 7- 
rays from embedded SNRs are produced by 7r° decay would be extremely 
precious to confirm the validity of the entire scenario. Such an evidence may 
come from the detection of neutrinos coming from the same sources. 

Several Neutrino telescopes are already operating, or under construe- 



tion, around the world (see e.g. [IHl)- Due to their large detection vol- 
ume and observation time, these instruments may have the capability to 
detect at least few neutrinos from some SNRs, especially if they are embed- 
ded in a dense medium. It is important to observe that the energy thresh- 
old above which Neutrino Telescopes may identify astrophysical neutrino 
sources ( > 100 GeV) is approximatively the same above which atmospheric 
Cherenkov telescopes can detect photons. This coincidence offers the inter- 
esting possibility to directly test if HE photons and neutrinos have the same 
hadronic origin confirming the half-century pending hypothesis that HECRs 
are accelerated by SNRs. Having that purpose in mind, in this paper we 
estimate the expected neutrino/photon flux ratio in the TeV energy range 
from embedded SNRs. 

Our computation of the Z//7 ratio from pion decay is similar to that 
followed by the authors of ^T] for RXJ1713-3946 (see also ^21) and confirm 
their main results, a part small differences which we discuss in SecOl 

In SecJUwe discuss the production of TeV antineutrinos by neutron decay 
in GMC complexes. In jT^] it was proposed that this process may give 
rise to a detectable antineutrino flux from the Cyg 0B2 association and 
from the GC That analysis was based on Ultra High Energy CR anisotropy 
measurements which, however, have a very low statistical significance or, in 
the case of the GC, have not been confirmed by other experiments. We will 
show that unless the primary particle spectrum is considerably different from 
that expected from ordinary shock acceleration, an observable contribution 
of neutron decay to the total neutrino flux can hardly be expected. Finally, 
in SecElwe apply the results of SecElto estimate the muon neutrino flux from 
the MC complex in the GC on the basis of the new HESS measurements [Zj. 

2 Photon emissivity by ttq decay 

The process under consideration is pp -^ NN + Htj-it'^''^ where n.^ is the pion 
multiplicity and N is either a proton or a neutron. The general expression 
of the photon emissivity is 
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where we assumed that hydrogen, having density nn, is the main proton 
target. When the proton energy is much larger than the proton mass, i.e. 
Ep > several GeV, the differential pion production cross section is well ap- 



proximated by a scaling expression 
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where x = E^o/Ep and ctq — 3 x 10~ cm~ . Using the numerical package 
PYTHIA pS] we verified that the expression of the scaling function adopted 
by the authors of il4 , fno{x) = 0.67(1 — x)^'^ + O.Se"^^"^, provides a good fit 
of simulated pp scattering data up to 500 TeV. 

In all cases in which the proton spectrum is well approximated by a power 
spectrum 

np{Ep) = no(^j , (3) 

the vr" emissivity becomes 

Q^oiE^o) ~ np(E°)ao c nnYia) , (4) 

where 

Y{a) = / dxx°'~'^fT,o{x,a) (5) 

Jo 

Hence the photon emissivity is 
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where we used E'^''{E^) =E^(l + 



It is also useful to write the mean value of the ratio between photon 
energy and that of the parent proton {photon elasticity), which is 

1 r [^ 

(E^/Ep) = 77 / dcos9 (1 — cos6') / (ia;/^o(a;) 

2 J -TV Jo /J\ 



/ dx f.,ro{x) = 0.18 
Jo 



'0 

where 6 is the angle of the photon momentum in the lab frame respect to 
that of the primary proton. We conclude this section by observing that by 
writing the proton spectrum in the form Eq. (jH)), we implicitly assumed that 
the ultra-violet cut-off is well above the maximal energy relevant for 7-ray 
observations. This is well motivated by the fact that the energy spectra 
observed by HESS for RXJ1713. 7-3946 and the GC jTj, which we assume 
to be representative of all embedded SNRs, are steady power laws up to 
~ 10 TeV. 



3 Muon and Electron Neutrinos from charged 
pion decay 

Inelastic pp collisions lead to roughly the same number of vr'^'s, vr+'s and 
7r~'s . If , as we are assuming here, 7r° decay is the dominant source of high 
energy photons, a comparable emission of z/^ and Uf, has to be expected. The 
emissivity of muon neutrinos and antineutrinos produced by direct vr"*" and 
7r~ decay is 

m^ f°° rlW 

Ql(E.,) ^ -^^ / -^QME.) , (8) 

where Qt,±{Et,) = Q^o(_E^) x (1 ± e), e being the asymmetry between the 
positive and negative pions effectively produced in pp collisions: 

(9) 
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where r = 771^/777,7^. Since we don't know the details of the interaction, the 
parameter e cannot be evaluated analytically and we have to recourse to 
numerical simulations. We performed a simulation of proton-proton interac- 
tions up to an incident energy of about 500 TeV with pythia [15 . Our best 
fit for e is 

e(s) = 0.1 - 6 • 10"^ ln(s/so) - 0.15 ■ ^/s^s , (11) 

with s = 2r7ip + 2m,pEp and sq — 1 GeV^. In Fig Q we show the energy 
dependence of our e as well as our best fit. We can now write 

QliE^) - irh^y " ^-^Q.iE^) , (12) 

where (e) represents the average of e over the relevant energy range which, 
in our case, is 10 GeV — 500 TeV. The computation of the contribution of 
secondary muon decay /i -^ ev^Ve to the neutrino emissivity is slightly more 
involved. The contribution of this process to the neutrino emissivity is 

QUE.) = /'" dE, QUE,) ^EiE^ , (13) 
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Figure 1: Energy dependence and best fit of e, as provided by pythia. 



where from the kinematics it follows -E™'" — E^ and E^^^ can be safely posed 
to oo if the muon spectrum is broad enough. The muon emissivity is 
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The neutrino emissivity can thus be estimated as done in ^| by using the 
spectrum-weighted momenta 
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where (l/°~^)o and {y"'~^)i are given in the same textbook for the different 
flavours. The total neutrino emissivity by both positive and negative pion 
decay is given by the sum of (fT^ and of (|TH|l . 

We can account for kaons contribution as well, provided that we can 
relate their flux to that of pions. Again, a simulation performed with the 
help of PYTHIA allows us to determine that the mean ratio K/n is about few 
percent. Thus, accounting for the well known branching ratio for the decay 
K -^ fiu, that is 63.5%, ^ and keeping in mind that no signiflcant dynamical 

^We neglect here three-body decay which have only a few percent BR. 



differences arise between pions and kaons, we found the correction to tfie 
neutrino emissivity due to pion. Tliis come out to be few percent. 

Finally, using the expressions for the spectrum-weighted momenta and a 
mean asymmetry e of about 5% we find the values of the ratios QuJQ^ which 
we shown in Tab^for several values of the spectral index a. 
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Table 1: Ratio between neutrino and 7-ray emissivities for several values of 
the primary proton spectral index a. 

These values are close, though not coincident, with those given in [TT] . 
The percentage difference ranges between 10% and 30%, being maximal for 
electron antineutrinos. Most likely, such a discrepancy is due to a value 
of the pion charge asymmetry that in ^T] is larger than the value adopted 
here. We think that the value of e used in jTT] is correct only at low energies 
(~ 100 GeV) and that the authors may have not taken into account the 
decreasing behaviour of e at larger energy. Indeed, from Fig.l we see that the 
scaling behaviour is reached only above several hundred GeV. Nevertheless, 
it should be noted that being the uncertainty in the value of e provided 
by PYTHIA of the order of 20 %, the values of the neutrino-photon relative 
emissivities given in Tab. ^ are almost compatible with those given in ^I] 
within errors. 



3.1 The effect of Neutrino oscillation 

As well established by several experiments, neutrinos undergo flavor oscilla- 
tions during their propagation in vacuum. Since the typical SNRs distance 
from the Earth is much larger than the neutrino oscillation length around 
the TeV, the phase of oscillations is very large so that we can safely deal with 
averaged vacuum oscillations. In this case the flavor oscillation probability 
can be written 
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Table 2: Oscillation probability in the average vacuum oscillation hypothesis, 
calculated with the mixing angles reported in the text. 



where Viv is the unitary mixing matrix among the lepton flavours Z, /' 
e, /i, r. After the propagation the neutrino flux is 
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E^.^ 
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F^ being the expected neutrino flux in the absence of oscillations. 

Several experiments have been performed, and are now in progress, to 
determine the values of the mixing parameters. Adopting the standard 
parametrization ^Tj for the mixing matrix ?7jj, we can describe our present 
knowledge as: 

• 012 ~ 00 ~ (33.9+1^)° by solar and KamLAND data [TBI. 

• 6*23 ~ Qatm — 45° ± 10° by atmosferic data (see e.g. [H]). 

• 013 ~ ± 10° from CHOOZ data [IS]. 

The CP violating phase dcp is still unknown and we assume here that it 
is negligible. Under this hypothesis the oscillation probability matrix Pij is 
symmetric and equal for both neutrinos and antineutrinos. We show it in 
TabEI 

Using Eq. (fT7|l we calculate the Z//7 ratios after propagation which we 
shown in Tab. 3 As an effect of oscillations the flavor composition is almost 
isotropized: the expected fluxes for each flavor are the same within a 20%. 
Another important effect is the appearing of a z/^ component, even if it is not 
generated at the source. 

Neutrino telescope looking for up-going muons produced in the Earth are 
mainly sensitive to the incoming flux of v^ and v^ . As a reference we write 
here the relation between the arrival v^^ + i/^ and 7-ray fluxes from a SNR 
accelerating primary protons with a power-law spectrum with a = 2.2 



XK 



|a=2.2 



0A6F^{E^ 



We will not discuss here the errors introduced by the uncertainties in our 
knowledge of the mixing angles. As showed in jirj, they are much smaller 
than the expected experimental errors. 
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Table 3: Z//7 ratios as seen at Earth after propagation. 

It is important to notice that since "averaged" oscillations are energy 
independent they do not affect the power-law energy spectrum expected for 
neutrino at the source. As it has been pointed out in 20^ , a spectral distortion 
in the flux of astrophysical neutrinos might be the signal of long-wavelength 
energy- dependent oscillations, possibly due to the existence of new families 
of sterile neutrinos. 



4 Antineutrinos from neutron decay 

The interest about this antineutrino production mechanism has been stim- 
ulated by claims of possible detections of CR anisotropy in the EeV range. 
Namely, the HIRES f2l1 and AGASA [22] collaborations found an excess of 
CRs with energy close to the EeV coming from the direction of Cygnus-0B2, 
though with low statistical significance. CR excesses was also detected by the 
AGASA |22] and SUGAR J2H1 experiments in directions close to the Galactic 
Centre (GC) and in the same energy range. Even if the existence of this 
anisotropy has not been confirmed by AUGER |21], theoretical arguments 
suggest that a positive signal may be found when the statistics of that experi- 
ment will be improved J2S|- Both the Cyg-0B2 star association, laying in the 
Galactic plane at a distance of 1.7 kpc from the Earth, and the region within 
100 pc from the GC, at a distance of about 8.5 kpc, are known to be active 
star forming regions rich of molecular gas and dust. The origin of the CR 
anisotropy in their directions have been imputed to a neutron emission from 
those systems with energy up to few EeV. At that energy, in fact, neutrons 
can travel over galactic distances and, differently from protons and compos- 
ite nuclei which are isotropised by the galactic magnetic field, they should 
keep the information about the angular position of their source unspoiled. It 
was proposed that if primary nuclei are accelerated above the EeV in that 
system, their PD onto the Far Infra Red (FIR) radiation which is present 
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in Cyg-0B2 and in the GC regions could efficiently produce neutrons of the 
required energy ^3]. The presence of a FIR background is a common feature 
of star forming regions. It is due to the absorption and reprocessing of optical 
and UV stellar radiation by the thick dust which is generally present in those 
systems. The spectrum of this radiation is almost that of a black-body with 
temperatures typically ranging between 10 and 100 K. 

The authors of jT3| suggested that, would the detection of a neutron flux 
from Cyg-0B2 and the GC be con&med, a "guaranteed" flux of antineutrino 
in the TeV energy range has also to be expected from those sources. Their 
argument is based on the hypothesis that the primary nuclei spectrum be- 
tween 1 PeV and few EeV is a rather steep power law (they assumed a > 3 
as for the CR spectrum observed on the Earth) and that TeV antineutrinos 
are produced by the decay of neutrons, with energy close to the PeV, which 
are produced by the PD of nuclei onto the UV radiation in the proximity of 
massive stars. Although this is an interesting possibility, it should be noted 
that the neutron spectrum was not measured and that, since neutrons should 
be produced in the nearby of the source and do not diffuse in the Galactic 
magnetic field, it is questionable that its slope coincides with that of ordi- 
nary CRs. Diffusion of primary nuclei in the MC region should be almost 
absent since, even assuming a magnetic field strength of ~ 100 fiG, their 
giration radius at energies above the EeV is comparable to the cloud size. 
Furthermore, only few, if any, SNRs may be able to accelerate nuclei beyond 
the EeV. 

Here we consider a more general scenario in which neutrons may not be 
produced with the energy needed to reach the Earth and be detectable and 
try to relate the antineutrino flux from neutron decay to the 7-ray flux which 
may be observed coming from the same source. We will determine the neu- 
tron production rate in the nearby of the SNR by applying some of the results 
derived in a previous work ^^. As secondary neutrons can be produced not 
only by the PD of high energy nuclei but also by hadronic collisions, mainly 
P + Pgas -^ n + p + pions, both processes need to be considered. We will show 
that, although PD is expected to be the most effective production process 
of EeV neutrons in MCs, at lower energies the hadronic scattering channel 
should be the dominant neutron production channel. 

4.1 Neutrons from nuclei photo-disintegration 

SNRs are expected to accelerate nuclei by a ffist order shock acceleration 
mechanism. In the linear shock acceleration theory all nuclear species are 
expected to have approximately the same power-law energy spectrum up to 
some UV cut-off depending on the charge number Z. 
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Above a centre of mass energy of a few MeV, a nucleus can photo- 
disintegrate into a lighter nucleus and one or more nucleons, due to the 
resonant giant dipole interaction with background photons. This process 
was studied in details in [SEl I23- Starting from the results presented in 
those papers, some of us derived the neutron production rate due to the PD 
of nuclei onto a thermal background of photons at the temperature T j23] . 
This is 

T^E, T) ^ 4 X 10-- U^AE, T) (j^^) s"^ , (19) 

where ^a is an order one parameter depending on the nuclear species. The 
function $^(i5,T), which has been defined in j23], depends on the product 
ET only. It is peaked at the energy 

where it takes a value depending on the nuclear species. For Iron nuclei 
(^^Fe) $56(-E'*) ~ 1 , it is an order of magnitude smaller for ^^O and ^^C, and 
3 X 10^^ for ^He. Although '^He are the most abundant composite nuclei, 
heavier species, especially Iron, give a comparable contribution to the PeV 
neutron production yield due to their larger PD cross section. 

Equation ()19|) has been derived under the assumption that the photon 
distribution is that of a black-body, or a grey-body (diluted black-body), 
with temperature T. Since each neutron takes, on average, a, 1/A fraction 
of the parent nucleus energy, it follows from pn|) that PeV neutrons will be 
mainly produced by nuclei PD onto a photon background with T ~ 1 eV. 
However, due to the huge opacity of the dense molecular gas practically no 
optical and UV photons should be present in GMC but in small H II blisters 
produced by the photoionization of the hydrogen gas around massive stars. 
The radius of these regions, which are generally known as Stromgren spheres, 
is expected to be roughly (see e.g. |2H]) 

R ^02 ( ^ \ '" ( ^ \ '"' ( ^^ ^ '"" 

Uhu-U. |^^q48s-i^ V2-6 X 10-23 cm-3s-V Vl03cm-3y P"^' 

(21) 
where F is the UV photon production rate of OB stars and /3 is the re- 
combination coefficient. The size of the Stromgren spheres grows relatively 
little during the short main sequence life-time of these massive stars |28j . 
Furthermore, the contribution of dust to the opacity may further reduce 
i?Hii- In dense GMCs, having uh = 10^ — 10^ cm"^ a characteristic size 
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of 10 pc and a population of hundreds of OB stars, the H II phase fill- 
ing factor is smaller than 10~^. Since the photon density in each blister 

J? 
is TT-uv — — 52 — — ^^^ cm"^, the overall mean UV photon density in a 

dense MC is roughly n^y ~ 1 cm~^. An order of magnitude larger photon 
density may be found in rich systems like the Cyg-0B2 association. As a 
consequence, in such a kind of systems, the H II phase filling factor is close 
to unity and the UV photon density nuv ~ 10 cm~^ which is in agreement 
with the value adopted in [T3]. Using ® and (jT^ we find 

Qn{E^) = J2fA np{En) Ta{AE^,T) 

(22) 

where Ja is the relative abundance of all the nuclear species with atomic 
weight A in the SNR which undergo effective PD. For the proton spectrum 
we used the unit GeV~^cm~^. 

4.2 Neutrons from p-p scattering 

Besides the PD process discussed in the previous section, EeV neutrons might 
also be produced by the collision of UHE nuclei with the dense gas in the 
MC environment. The relevant neutron production channel is the charge 
exchange pp inelastic collision pp^^s —>■ pn + "oIj^tc, where m^^ is the pion 
multiplicity. In the case of a power-law spectrum for the HE protons the 
neutron emissivity due to this process was found to be |E| 

Ql^iEn) = nHnp{En)ao m^ c K I{a) 

i5„/i:.A / ^H \fl{a)\ ^,,-i_-3„-i (23) 



- 2 X 10- «,(£„) ^^jj^j ^^ j GeV-'c,„ 

where K = 4/5 is the inelasticity, and m„ ~ 0.24 the neutron multiplicity. 
We neglected here a weak logarithmic energy dependence of the cross section. 
The function I (a) is defined by 

I{a) = f dxx''-^h{x) (24) 

Jo 

where h{x) ~ 0.064 (1 - xf + 0.094x^1 - x and x = E^/Ep. 
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By comparing ()23|) with (j22j) it is evident that under most common condi- 
tions in GMCs, pp scattering should be the main neutron production channel 
at energies around the PeV. 

An approximated expression for the antineutrino emissivity from neutron 
decay can then be derived by treating the neutron decay factor as a step 

((1771 \ 
— En ] and replacing the antinuetrino energy distribution 
Tn J 

by a delta-function peaked at the mean value Ejj ~ 6 x 10 ^ rrin (here we 
followed the same approach adopted in [T3] )• By doing that, we find 

g^(E^,x) = 0.1 ^ (^) " g,(i?^,x) . (25) 

Comparing this result with that given in SecEl it is clear that, unless the pri- 
mary proton spectrum has a quite unusual behaviour, neutron decay should 
give only a negligible contribution to the total neutrino emissivity of an em- 
bedded SNR. 

The relative contribution of neutron decay might somewhat be enhanced 
if the proton energy distribution in the MC has a strong spatial dependence. 
This may be indeed the case for protons which are shock accelerated by a 
steady SNR, as the expected energy spectrum upstream the shock front is 

m 

n,{E,r) = n,{E,0) exp {"^^j • (26) 

Here r is the distance from the SNR shock; no(-E) is the proton spectrum 
downstream the shock, which we assume to be given by Q ; us is the shock 
velocity; D{E) is the diffusion coefficient; k is an order one coefficient. For 
the sake of simplicity we also assumed that us and D{E) do not depend 
on r. By expressing the energy dependence of the diffusion coefficient in 
the usual form D{E) = D{Eq){E/Eo)^, the meaning of (jSHI) become quite 
evident: upstream the shock the proton spectrum has an exponential low- 
energy cutoff at 

^„,.(,) = E.(^) . (27) 

Assuming Bohm diffusion in the acceleration region (note that adopting a 
larger diffusion coefficient in that region would imply a too low maximal 

acceleration energy), i.e. D{E) = - c rL{E), where riiE) is the Larmor 
radius, we find 

-'^) == "> (riv) (500^) "' (lo^) " "- <^«' 
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It is evident that PeV protons can diffuse into the cloud much more deeply 
than 10 TeV, so that secondary neutrons, as well as TeV 77^ from their decay, 
are produced in a volume larger than that in which photons and neutrinos 
are emitted by hadronic scattering. Note, however, that the characteristic 
distance over which protons propagate before loosing their energy due to pp 



scattering (rioss ^ (c uhcto) ^), is 



1/2 



(29) 
Therefore, energy losses significantly reduce the effective propagation volume 
of primary protons in those cases in which the gas density and magnetic field 
strength are very large. Unfortunately, these are also the most interesting 
cases from the observational point of view. At the end, we found that in all 
relevant physical situations the volume gain factor is not sufficient to make 
the enhancement of the contribution of neutron decay to the total neutrino 
emissivity observationally relevant. 

We conclude this section by stressing that all the previous considerations 
apply only if the particle primary spectrum is a power law. This will not be 
the case, for example, if primary nuclei are accelerated by a strong voltage 
drop around a fast rotating neutron star. In that case, in fact, primary nuclei 
may have an energy spectrum peaked at the PeV, or above, so that their PD 
onto stellar radiation be dominant respect to hadronic scattering and neutron 
decay be the main neutrino production channel. 



5 Neutrinos from the Galactic Centre 

In this section we apply our previous results to estimate the neutrino flux 
from the dense molecular cloud complex in the Galactic Centre region. An 
extended very-high-energy 7-ray emission has recently been discovered by the 
HESS coming from a region roughly delimited by |/| < 0.8° and |6| < 0.3° 
in galactic coordinates [Zj. The coincidence between the angular position of 
the centroid of this emission and that of the dynamical centre of the Galaxy 
suggests that both are located at the same distance from the Earth, which 
is about ~ 8.5 kpc. The energy spectrum of this source was measured to 
be (1.73 ±0.13 ±0.35) x ^~(2.29±o.o7±o.o2) TeV-^cm-^s-^sr-^ Superimposed 
to that diffused emission HESS found a point-like source (J1745-290), which 
was also early detected by HESS [SlI (and confirmed by MAGIC |31j) with 
an energy spectrum (2.50 ± 0.2) x Er^^^ ' TeV~Vm^^s^^, in a position 
compatible either with the supermassive black hole in Sgr A* or with the SNR 
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Sgr A East. This source was subtracted when determining the spectrum of 
the extended emission. 

As the authors of ^ suggest, the detected 7-rays are hkely to be produced 
by the hadronic interaction of locally produced cosmic rays with the dense 
molecular gas. A relevant support in favour of this hypothesis is given by the 
strong correlation found between the 7-ray and the millimetric emissions from 
the CS (Carbon Sulfide) in that region, that compound being a well known 
tracer of molecular hydrogen. The coincidence, within errors, between the 
spectral slope of the GC ridge emission detected by HESS and that of J1745- 
290, suggests that the primary particles responsible for both signals may 
have the same origin. Indeed, the total CR energy required to explain both 
of them is compatible with that typically expected from shock acceleration 
by the remnant of a single SN. These considerations point to Sgr A East as 
the most plausible primary source (see jHSl 12^1 about the identification of 
J1745-290 with that SNR). 

This elegant scenario may find further support would neutrinos from the 
same region be detected. The expected neutrino flux can be easily deter- 
mined by applying the results of Sec. El It has to be noted that due the 
relatively poor angular resolution (of order < 1°) achievable with Neutrino 
Telescopes and the very low expected neutrino fluxes, it will be quite hard, 
if not impossible, to disentangle the neutrino point-like emission from that 
coming from the Galactic Centre ridge. Thus, we have to sum the neutrino 
flux from both sources. This is straightforward since, as the 7-ray spectral 
slopes of these sources coincide within errors, the same will be true for neu- 
trinos. We use here the central value of the slope measured by HESS for the 
GC ridge, i.e. a = 2.29. The total 7-ray flux from the GC ridge source is 
determined by multiplying the diffused flux measured by HESS by its solid 
angle opening, i.e. 4 x 0.8° x 0.3° = 3 x lO^'* sr. It is about the double of the 
flux from J1745-290. Therefore, the total 7-ray flux, above the TeV, coming 
from the GC region is roughly (we only write central values in the following) 

(E \ '^'^ 
p^ j TeV-^ cm-' s'' . (30) 

Using our results of SeclHl we derive the Z//7 ratios for a = 2.29 at the source 
and at the Earth are shown in TablUfor the different neutrino species. We 
notice here that the TeV 7-ray attenuation due to pair-production scattering 
onto the IR photon background should be negligible. In fact, whenever the 
radiation temperature is below 100 K, the mean IR photon energy implies a 
7-ray threshold energy well above 10 TeV. Gamma-ray scattering onto the 
UV and X-ray backgrounds present in that region was also showed to give 
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a = 2.29 


at source 


at Earth 


^e 


0.24 


0.23 


Ve 


0.21 


0.21 


^H 


0.43 


0.22 


v^. 


0.43 


0.22 


Vr 





0.22 


Vr 





0.22 



Table 4: Z//7 ratios before and after oscillations. 

rise to a negligible attenuation [H!^ . 

It is now straightforward to derive the total expected niuon neutrino flux 
which is 



i^GC 



XEu 



3.3 X 10 



-12 



3l. 

TeV 



-2.3 



TeV ^ cm ^ s 



(31) 



This flux is about the double of that estimated in 32^ using only the 2004 
HESS data [31j. Since the detection of secondary up-going muons produced 
by neutrino interaction in the Earth is the most promising strategy to identify 
TeV neutrino from astrophysical sources, Neutrino Telescopes in the North 
hemisphere are best suited to look for a signal from the GC as it is located 
in the southern Sagittarius constellation. Using the v^ + z7^ flux given in 
(jnij, and the ANTARES effective area for up-going i/^s given in [33], we 



estimated the expected detection rate by this experiment to be 0.07 yr~^ 
which is, unfortunately, quite small. A km Neutrino Telescope to be built 
in the Mediterranean sea (e.g. NEMO [221 o'^ NESTOR [2Zj upgrading) has 
better chance to get a positive detection. Indeed, since the effective area of 
such an instrument should be at least 20 times better than that of ANTARES, 
the expected rate is 1.5 yr~^. As the expected background rate of such device 
is ~ 0.3 ev/yr . [32 , the observation time required to detect that source at 
95%C.L is about three years. 

6 Conclusions 

In this paper we estimated the neutrino/photon flux ratio in the TeV energy 
range from embedded SNRs. We did that under the hypothesis that SNRs 
shock accelerate protons and composite nuclei with a power law energy spec- 
trum at least up to the PeV. By assuming that all 7-rays of that energy are 
produced by vr" decay, we flrst computed the relative emissivities of muon 
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and electron neutrinos produced by charged pion decay. Then, we took into 
account for vacuum oscillations to determine the 1/4/7 Aux ratio to be ex- 
pected on the Earth for all lepton families. Our results are similar, though 
not coincident, to those previously derived by other authors [HIISZI- We also 
considered the possible contribution that neutron decay may give to the TeV 
antineutrino flux from a GMC embedding an active SNR. Using the results 
which we derived in a previous work [23], we showed that in dense molecular 
complexes (e.g. in the GC) hadronic scattering should be the dominant neu- 
tron production mechanism at the relevant energies (~ PeV). That channel 
was not considered in [13]. Unfortunately, we found that, although such a 
neutrino production channel is less severely constrained by 7-ray observa- 
tions, it is quite unlikely that under reasonable assumptions it could add a 
significant contribution to the total neutrino flux from such kind of sources. 
As an application of our results we used the recently released data of 
the GC observation campaign by the HESS collaboration [2] to estimate the 
muon neutrino flux from that region. We found that the km^ Neutrino Tele- 
scope to be built in Mediterranean sea may have some chances to detect 
that source. Although the expected neutrino detection rate is rather small, 
the possible observation of a positional coincidence between the 7-ray and 
the neutrino source as well as the measurement of a flux ratio at the level 
estimated in this paper, would add significant evidence in favour of the hy- 
pothesis that hadron acceleration is taking place in that region. Furthermore, 
since embedded SNRs should be quite common objects in the Galaxy and 
VHE 7-ray and neutrino astronomy have just started, other sources of that 
kind may be discovered in both channels in a not too far future. 

Acknowledgements 

We thank V. Flaminio, T. Montaruli, T. Sjostrand and F. Vissani for several 
valuable discussions. The authors are members of the ANTARES collabora- 
tions. L.M. thanks G. Usai for help using the Pythia package. 

References 

[1] V.L. Ginzburg and S.I. Syrovatskii, Origin of the Cosmic Rays, Pergamon 
Press, London, 1964. 

[2] T. Montmerle, Astrophys. J. 231 (1979) 95. 

[3] F.A. Aharonian and A.M. Atoyan, Astron. Astrophys. 309 (1996) 917. 



17 



[4] S.D. Hunter et al, Astrophys. J. 481 (1997) 205. 

[5] D. F. Torres, G. E. Romero, T. M. Dame, J . A. Combi and Y. M. Butt, 
Phys. Kept. 382 (2003) 303 arXiv:astro-^/0209565j . 

[6] F. Aharonian et al. [The HESS Collaboration], Nature, 432 (2004) 75. 

[7] F. A. Aharonia n [HESS Collaboration], Nature 439 (2006) 695, 
|arXiv:astro-ph70 603021 . 



C. Heiles and R. Crutcher, arXiv:astro-ph/0501550[ 



[9] F. Yusef-Zadeh et al, Astrophys. J. 466 (1996) L25 . 



[10 
[11 

[12 

[13 
[14 
[15 
[16 
[17 

[is; 

[19 
[20 

[21 



T. Montaruh, Eur. Phys. J. A 24S1 (2005) 103. 

M. L. Costantini and F. Vissani, Astropart. Phys. 23 (2005) 477 
[arXiv:astro-ph/0411761. . 



M. L. C ostantini and F. Vissani, AIP Conf. Proc. 794 (2005) 219 
|arXiv:astro-ph/0508152, . 



L. A. Anchordoqui, H. Goldberg, F. Halzen and T. J. Weiler, Phys. Lett. 
B 593 (2004) 42 |arXiv:astro -ph/0311002| . 

P. Blasi and S. C olafrancesco, Astropart. Phys. 122 (1999) 169 
|arXiv:astro-ph/9905122| . 



T. Sjostrand, L. Lonnblad, S. Mrenna and P. Skands, 



|arXiv:hep-ph/0308153 



T.K. Gaisser, Cosmic Rays and Particle Physics, Cambridge University 
Press, 1990. 

S. Eidelman et al. [Particle Data Group], Phys. Lett. B 592 (2004) 1. 

B. Aharmim et al. [SNO Collaboration], Phys. Rev. C 72, 055502 (2005) 
|arXiv:nucl-ex/0502021| . 



M. Apollonio et al, Eur. Phys. J. C 27, 331 (2003) 
|arXiv:hep-ex/0301017| . 



R. M. Cr ocker, F. Melia and R. R . Volkas, Astrophys. J. Suppl. 141 
(2002) 147 |arXiv:astro-ph/0 106090 . 

D. J. Bird et al. [HIRES Collaboration] , arXiv:astro-ph/9806096J 

18 



[22 
[23 
[24 
[25 
[26 

[27; 

[28 

[29 
[30 

[31 

[32 

[33 

[34; 

[35 
[36 



N. Hayashida et al. [AGASA Collaboration], Astropart. Phys. 10 (1999) 
303 arXiv:astro-ph/9807045 1 . 

J. A. Bellido, R. W. Clay, B. R. Dawson and M. Johnston-Hollitt, As- 
tropart. Phys. 15 (2001) 167 |.arXiv:astro-ph/0009039j . 

A. Letessier-Selvon [Pierre Auger Collaboration], 29th International 
Cosmic Ray Conference, Pune, 2005, |arXiv:astro-ph/0507331, . 

P. Grasso and L. M accione, Astropart. Phys. 24, (2005) 273 
|arXiv:astro-ph/0504323| . 



J. L. Puget, F. W. Stecker and J. H. Bredekamp, Astrophys. J. 205 
(1976) 638. 

F. W. Stecker and M H. Salamon, Astrophys. J. 512 (1992) 521 
|arXiv:astro-ph/9808110| . 



R.I. Diaz-Miller, J. Franco and S.N. Shore, Astrophys. J. 436 (1994) 
795. 

M.A. Malkov and L. O'C. Drury, Rep. Prog. Phys. 64, (2001) 429 

F. Aharonian et al. [The HESS Collaboration], Astron. Astrophys. 425 
(2004) L13 |arXiv:astro-ph/0408145| . 

J. Albert et al. [M AGIC Collaboration], Astrophys. J. 638 (2006) LlOl 
|arXiv:astro-ih70512469| . 



R. M. Crocker, F. Meha and R. R. Volkas, Astrophys. J. 622 (2005) L37 
|arXiv:astro-ph/041147~ . 



R. M. Crocker, M. Fatuzzo, R. Jokipii, F. Melia and R. R. Volkas, 
Astrophys. J. 622 (2005) 892 |arXiv:astro-ph70 408183 . 

F. Aharonian et al. [The HESS Collaboration], Astron. Astrophys. 425 
(2004) L13 |arXiv:astro-ph/0408145| . 



E. Aslanides et al. [ANTARES Collaboration] , |arXiv:astro-ph/9907432 



C.N. De Marzo [NEMO Collaboration], Nucl, Phys. Proc. Suppl. BlOO 
(2001) 344. 



[37] S. E. Tzamarias [NESTOR Collaboration], Nucl. lustrum. Meth. A 502 
(2003) 150. 



19 



